Frequency Dependent Specific Heat from Thermal Effusion in Spherical Geometry. 
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We present a novel method of measuring the frequency dependent specific heat at the glass 
transition apphed to 5-polyphenyl-4-ether. The method employs thermal waves effusing radially 
out from the surface of a spherical thermistor that acts as both a heat generator and thermometer. 
It is a merit of the method compared to planar effusion methods that the influence of the mechanical 
boundary conditions are analytically known. This implies that it is the longitudinal rather than 
the isobaric specific heat that is measured. As another merit the thermal conductivity and specific 
heat can be found independently. The method has highest sensitivity at a frequency where the 
thermal diffusion length is comparable to the radius of the heat generator. This limits in practise 
the frequency range to 2-3 decades. An account of the 3a;-technique used including higher order 
terms in the temperature dependency of the thermistor and in the power generated is furthermore 
given. 

PACS numbers: 64.70.P- 
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I. INTRODUCTION 

A hallmark of the glass transition of supercooled liquids is 
the time dependence of physical properties — also called 
relaxation. This time dependence is typically seen on 
the timescale of the so-called Maxwell relaxation time 
tm- tm is defined as 



TM 



Goo 



(1) 



where 770 is the low frequency limiting (DC) shear vis- 
cosity and Goo the high frequency limiting elastic shear 
modulus. The relaxation is often most easily studied in 
the frequency domain. At high frequencies the liquid 
shows solid-like elastic behaviour whereas at low frequen- 
cies it shows liquid-like viscous behaviour. In general the 
viscoelastic stress response to oscillatoric shear strain is 
described by the complex dynamic shear modulus G [J. 
The structural relaxation is also observed in the specific 
heat. This has long ago been studied in the time domain 
as enthalpy relaxation but in 1985 frequency domain 
studies or AC-calorimetry 01 were introduced in the 
research field of supercooled liquids. Common to these 
methods is that they excite thermal oscillations in the 
sample at a cyclic frequency to. For a heat diffusivity, 
D, this frequency is associated with a characteristic heat 
diffusion length, 



(2) 



the magnitude of which gives the range of the excited 
thermal waves. When comparing \Id\ to a given sam- 
ple size, L, one can classify AC-calorimetry experiments 
as to whether the sample is thermally thin (L ^ \Id\) 
or thermally thick (L ^ A typical value of D 

is O.lmm^/s giving \Id\ « 0.1mm for a frequency of 
1/ = uj/{2tt) — IHz. This means that thermally thin 



methods are restricted to rather low frequencies [1, Q or 
very thin samples [5, 6]. In these methods the temper- 
ature field is homogeneous throughout the sample and 
the frequency dependent specific heat is derived directly. 
The AC-technique allows for easy corrections of the in- 
fluence of heat leaks in a non-adiabatic configuration. 
AC-calorimetry on thermally thick samples on the other 
hand has the advantage of covering a large frequency span 
[1, 01 • In this case the temperature field is inhomogenous 
and what is really measured is the effusivity, e = •\/Ac, 
from which the specific heat is derived. The effusivity 
expresses the ability of the liquid to take up heat from 
parts of its surface and transport it away (thermal effu- 
sion). This ability is dependent on two constitutive prop- 
erties, the heat conductivity. A, and the specific heat (pr. 
volume) , c. The specific heat can only be found to within 
a proportionality constant in effusion experiments from 
an (infinite) plane. However taking boundary effects into 
account for finite size plane heaters 0, Q one may in 
principle find A and c independently. The intermediate 
situation where the heat diffusion length is comparable to 
the sample thickness has been exploited in the so-called 
two-channel ac calorimeter [l^, HH. This fine method 
aims directly at an independent determination of A and 
c. 

The specific heat entering the effusivity is widely taken 
as the isobaric s pec ific heat, Cp. However recently it has 
been shown |l^Tl4| that in planar and spherical geome- 
tries it is rather the longitudinal specific heat, c;, that 
enters. When Ms and Mt are the adiabatic and isother- 
mal longitudinal moduli respectively, ci is related to the 
isochoric specific heat, cy, by 



Ms 



(3) 



In contrast, when Ks and Kt are the adiabatic and 
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isothermal bulk moduli respectively, Cp is given by 

— ^V- 



(4) 



Now Ms = Ks + 4:/3G, and Mt = Kt + 4:/3G. This 
means that q will differ from Cp in the relaxation regime 
where the shear modulus, G, dynamically is different 
from zero, and Kg and Kt differ significantly. 

In this work we study the frequency dependent specific 
heat from thermal effusion in a spherical geometry. Spe- 
cific heat measurements in a spherical geometry have ear- 
lier been conducted in the thermal thin limit jH, but in 
this paper we explore the thermal thick limit. 

The justification is threefold 1) We want to realize an 
experimental geometry in which the thermomechanical 
influence on the thermal effusion can be analytically cal- 
culated 2) The spherical container of the liquid is 
a piezoelectric spherical shell by which it is possible to 
measure the adiabatic bulk modulus [l5| on the very same 
sample under identical conditions. 3) The heat conduc- 
tivity and specific heat can be found independently, as 
the introduction of a finite length scale allows for sepa- 
ration of the two. 

We also in this paper give an account of the 3uj detection 
technique of our variant. 



II. LINEAR RESPONSE THEORY APPLIED TO 
THERMAL EXPERIMENTS 



The experimental method used in this study is an effusion 
method. In such methods a harmonic varying heat cur- 
rent, Re(Pe'"*), with angular frequency uj and complex 
amplitude P, is produced at a surface in contact with the 
liquid; this heat then "effuses" into the liquid. At the sur- 
face a corresponding harmonic temperature oscillation, 
Re((5Te*'^*), with complex amplitude ST, is created. The 
current and temperature can be thought of as stimuli 
and response (but which is which can be interchanged at 
will), and the two are linearly related if the amplitudes, 
|P| and \6T\, are small enough. Hence it is convenient 
to introduce the complex frequency-dependent thermal 
impedance of the liquid 



ST 



(5) 



this impedance depends on the properties of the liquid. 

The terminology thermal impedance derives naturally 
from the close analogy between the flow of thermal heat 
and the flow of electricity [l^ . Temperature corresponds 
to electric potential and heat current to electrical cur- 
rent. The analogy also applies to the heat conductivity 
which corresponds to the speciflc electrical conductivity, 
and the specific heat which corresponds to electrical ca- 
pacity. 



As an example the thermal impedance of a thermally thin 
sample is given as 



•2^1iq,thin — 
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iujcV 



(6) 



where V is the volume of the sample. Again an equa- 
tion completely equivalent to the electrical correspon- 
dence between capacitance and impedance. 

In the classical case of a planar heater in a thermal thick 
situation the thermal impedance is found to be [3] 



1 



liq, planar 



AVicJcX 



(7) 



where A is the plate area. The lateral dimension W has 
to be large {W ^ |^£)(a;)|) for this formula to be exact 
and unfortunately correction terms that take boundary 
effects into account decays only slowly as Q. 



III. THERMAL EFFUSION IN SPHERICAL 
GEOMETRY 



Thermal effusion is possible in a spherical geometry as 
well as in a planar geometry. Imagine a liquid confined 
between the to radii, ri and r2. The imposed oscillat- 
ing outward heat current, Rc (P(ri)e'"^*) , at radius ri 
results in a temperature oscillation. Re ((5T(ri)e*"*) , at 
the same radius. The heat wave propagating out into the 
liquid gives rise to a coupled strain wave due to thermal 
expansion. The heat diffusion is thus dependent on the 
mechanical boundary conditions and mechanical proper- 
ties of the liquid. It is widely assumed that the transport 
of heat in a continuum is described by the heat diffusion 
equation 



iujST = DV^ST. 



(8) 



By the same token one might take the heat diffusion con- 
stant as 13 = X/cy if the surfaces at ri and r2 are clamped 
and D = X/cp if one of the surfaces is free to move. This 
common belief is not entirely correct. If shear modulus 
is nonvanishing compared to bulk modulus the descrip- 
tion of heat diffusion cannot be decoupled to a bare heat 
diffusion equation (T2| - fT3 . [l7| . Such a situation arises 
dynamically at frequencies where uj w ^/tm- The fully 
coupled thermomechanical equations have been discussed 
in detail, by some of us, in planar geometry [l^ and in 
spherical geometry In general the solutions are com- 
plicated but in the case of a spherical geometry where the 
outer radius r2 is much larger than the thermal diffusion 
length, <C r2 the solution becomes simple. In 

terms of the liquid thermal impedance, Znq, at the inner 
radius, ri, one has 



Zliq 



ST in) 



(9) 
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It is very interesting that in this thermally thick limit, 
the thermal impedance is independent of the mechani- 
cal boundary conditions [l^[lJ|- I-e. in the two cases of 
clamped or free boundary conditions mentioned before it 
is neither the isochoric nor the isobaric specific heat that 
enters the diffusion constant and the thermal impedance 
but the longitudinal specific heat. Since a closed spheri- 
cal surface has no boundary there are no correction terms 
like the ones to Eq. ([7]) discussed above for planar plate 
effusion. Notice from Eq. © that associated with a given 
radius, ri, is a characteristic heat diffusion time for the 
liquid 

Ti = ri'^ci/X. (10) 

c; and hence t; may in general be complex, but when c; 
is real, the imaginary part of Zuq peaks at w = I/t;. 



IV. Scje-DETECTION TECHNIQUE 

The SoJe detection technique we use (where the subscript 
e indicates that it is the electrical angular frequency not 
the thermal) is different from other 3cJe techniques (e.g. 
Ref. [3|) in two ways. Firstly, we include higher order 
terms in the temperature dependence of the thermistor 
used and in the power oscillations produced. Secondly, 
we use a direct measurement of the time dependent volt- 
age over a voltage divider instead of a lock-in amplifier. 



A. The principle. 

The Swc-detection technique makes it possible to measure 
the temperature with the same electrical resistor which 
generates the heat and thereby know the thermal current 
and the temperature at the same surface. The principle 
is in short the following: Applying an electrical current 
at cyclic frequency uic through a resistor a combined con- 
stant (zero frequency) and 2lOc component Joule heating 
is produced. This gives rise to a combined temperature 
DC-offset and temperature oscillation at 2a;o, the size of 
which depends on the thermal environment (the thermal 
impedance). The resistor is chosen to be temperature 
dependent, it is a so-called thermistor. The tempera- 
ture oscillations therefore creates a perturbation of the 
resistance at frequency 2uje- The voltage across the ther- 
mistor which by Ohm's law is the product of current and 
resistance thus contains a (and an additional lojc) 
component proportional to the temperature oscillations 
at 2ajc- In practise one rather has a voltage source than 
a current source. In order to get a 3^0 signal the ther- 
mistor is therefore placed in series with a temperature 
independent preresistor as detailed in the following. 

The result of the analysis is the frequency dependent 
thermal impedance at a frequency of 2aJo- It is conve- 
nient to express this thermal impedance as function of 



the frequency of the thermal signal which we designate 
Lj and notice that w = 2cjo. 

We first outline the simplest theory in order to elucidate 
the technique most clearly. However, to increase the sig- 
nal to noise ratio it is expedient to choose a thermistor 
with a large temperature dependence and to choose tem- 
perature amplitudes as high as possible although within 
the linear regime of the liquid thermal response. This ne- 
cessitates higher harmonics be taken into account. This 
detailed analysis is deferred to appendix [X] 



B. A detour on tlie complex notation. 

Complex numbers are known to be of great use in har- 
monic analysis of linear systems. Since the heat current 
is quadratic in the imposed electrical current and the 
whole idea of the 3a;o-detection technique is to exploit the 
inherent nonlinearity introduced by the temperature de- 
pendence of the thermistor the formalism becomes a little 
more troubled than in pure linear applications. Thus we 
have to drag around the complex conjugate terms in or- 
der to get product terms right. 

Generally a sum of harmonic terms, 

A = Aq + \Ai\cos{u!ct + (t>i) + ■ ■■ + \An\ cos{nujct + </)„), 

(11) 

with real amplitudes, Ao,|Afc|, and phases, (j}k, can be 
written as, 

A=^(^An+Ao + |Ai|e'("'=*+'^i) + \Ai\e-'^'^'^'+'*'''> + ... 
-f|^„|e*(""'=*+'^") + |yl„|e-'("'^'=*+'^"') . (12) 

By introducing the complex amplitudes Ak = \Ak\e^'^'' 
and the shorthand notation = e the sum becomes 

A = ^{Ao + AiEi + ...+ AnE„ + c.c), (13) 

where +c.c. means that the complex conjugated of all the 
terms within the parenthesis should be added, including 
the real constant term. Note that Eq = I, E^Ei = E^+i, 
and EkE* = Ek-i. 

Notice that any product of an n'th and m'th harmonic 
will produce both an n 4- m'th and n — mWi harmonic 
term. Furthermore the DC-components get doubled be- 
cause there is no difference between Eq and Eq. 



C. Fundamental equations for the voltage divider 

Consider the diagram of the electrical setup seen in Fig. 
[2 U{t) is the voltage of the source and V{t) the volt- 
age across the temperature independent preresistor, i?pre- 
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FIG. 1: Voltagedivider for the 3a;e-detection technique. The 
temperature-dependent resistor (the thermistor), R(T), is in 
thermal contact with the hquid probing its thermal impedance 
at 2uJe. The preresistor, Rpre, is temperature-independent. 



TABLE I: Experimental parameters. Top: Parameters char- 
acterizing the resistors in the setup, -Rprc is the resistance of 
the preresistor. Roc and Ta the characteristics of the ther- 
mistor according to Eq. (|19p ). Middle: Absolute size of the 
thermistor bead. Measured values indicates the eccentricity 
of the bead. The optimize value is the one use in the final 
treatment of the data, and is optimized to give best corre- 
spondence between the room temperature DC specific heat, 
and the literature value of this. Bottom: Applied amplitudes 
expressed as mean applied voltage, ((|C/i|)i/), and maximum 
temperature amplitude, (max(|r2 |)^,t). 



Rpic 


3884.30 


Roa 


0. 12510 


Ta 


2789.5K 


^1, Measured, long edge 


0.205mm 


"^1, Measured, short edge 


0.185mm 


'^l, Optimized 


0.2073mm 


mi). 


4.909V 


max(|r2|),,T 


2.253K 



U{t) and V{t) can be measured directly. R{T{t)) is the 
temperature dependent resistance of the heater/sensor 
(we will expUcitly state the temperature/time depen- 
dencies when appropriate, otherwise it is implicitly as- 
sumed) . 

The current /(<) through R{T{t)) and i?prc is given as 
and the power produced in the heater as 



Rim) 



P{t)^I{tfR{T{t)) = - 

{Rp,, + R{T{t))y 

Eq. (|14p gives the voltage divider equation, 



:U{tf. (15) 



R 



pre 



i?p,c + i?(T(t)) 



(16) 



which in principle allows for determination of T{t) from 
the measured voltages. 



D. The AC-method to lowest order 

The thermistor and the surrounding sample are placed 
in a cryostat that defines an overall reference tempera- 
ture, Tciyo. The power produces a (real) temperature 
change AT = T — T^yo at the heater/sensor depending 
on the thermal impedance, Z, (as defined in section |Tl]) 
the heater /sensor looks into. The temperature rise has 
a real DC-component, Tq, and complex AC-components. 
It can in general be written as 



AT = -{To+TiEi+T2E2+T3E3+T^E4 + ... + c.c.) (17) 



the dominant terms will be the DC component, Tq, and 
the second harmonic term, T2 (see next section). 

To first order in the temperature variations the thermis- 
tor has a resistance of 



R^ Ro{l + aiAT). 



(18) 



Here Rq and ai are considered constant at a given 
cryostate temperature but they vary with Tcryo- In fact 
the thermistor, which is based on a semiconducting ma- 
terial, has a resistance that over a larger temperature 
range is very well described by 



R = Ro 



(19) 



with Ta the activation temperature and R^o the infinity 
temperature limiting resistance, ai is related to Ta by 



1 



ai 



dR 



Ta 

cryo 



(20) 



i?(Tc,yo) dT 

Characteristic values of Rao and Ta are given in table |ll 
Substituting Eq. into Eq. we obtain 
1 _ 1 



V ^ 



l + ^(l + aiAT) 



U 



A+l 



{l-aAT)U (21) 



where A = Rg/Rp 



and the second expression 



is to first order in AT hke in Eq. HH]). 

In the simple case, which we first discuss, we assume that 
the function generator can deliver a pure single harmonic 
voltage, 



1 



U{t) = -(UiEi+c.c). 



(22) 



5 



Eqs. (ini) and (HH) can now be substituted into Eq. (^1)1 . 



V = (^l-a Q(To + T,E, + T2E2 + nE3+ 

T^Ei + ... + c.c.)]] liUiEi + U*El). (23) 



In principle the result contains all harmonics, but it turns 
out that the interesting ones are the first and third har- 
monic. The complex amplitudes are by inspection seen 
to be 



Vi 



1 



A + 1 



and 



^^3 = 



A + 1 



Ui - aToUi - -aT^Ul 



(24) 



(25) 



Notice that even if AT is written as an infinite sum (as in 
Eq. (flT)) ). Vi and V3 have only the included terms. This 
is a consequence of assuming a perfect voltage source. 

The T4 term can be neglected, since it is proportional 
to a higher order power term, P4, which is small. This 
approximation leads to the a use full relationship between 
the 3wc-term and the temperature amplitude at 2w, T2, 



2 (A + 1) 



T2U1. 



(26) 



Eq. dSll) and ^ can be solved for Tq and T2 if prior 
knowledge on the temperature dependency of the resis- 
tance of the thermistor exists, that is if A and a are 
known at the relevant temperatures. The T2 and Tq 
component terms are related to the power through the 
complex AC-thermal impedance, Z2 (the impedance at 
2we), and DC-thermal impedance, Zq as 



To = ZqPq and T2 ^ Z2P2- 



(27) 



The temperature dependency of the thermistor resistance 
can of course be found from calibration measurements. 
We have however found it to be more efficient to use 
an iterative solution technique eliminating the need for 
prior calibration. Eq. (p4| and ([26| are rewritten in the 
following form 



^ ^ [Ui - aTpUi - ^aT2U^) ^ 



T2 = -2 



aUi 



(28a) 
(28b) 



It can be seen that the equation for A consists of the 
O'th order expression A « ^^/Vi — 1 which would be valid 
if the temperature of the thermistor bead was the same 
as the cryostat temperature. This is not the case as the 
thermistor generates self heat, which is what the two ad- 
ditional terms correct for. 



A and T2 can be found by iteration in the following way. 
A first approximation of A is found by setting Tg = and 
T2 = in Eq. (|28ap . Based on this an initial estimate of 
ai is found from the temperature dependence of A. An 
estimate on T2 (and hence Z2) is afterwards found from 
Eq. (j28bp . By extrapolating Z2 to zero frequency an es- 
timate of To is found ^ (utilizing Tq = ZqPo). The 
estimated values for Tq,T2 and a can then be inserted 
into Eq. (|28ap giving a better estimate for A. This iter- 
ation procedure can be repeated a number of times and 
converges rapidly to a fixpoint for A and T2. 



E. Thermal-power terms 

If the 1. order approximation of the resistance of the ther- 
mistor, Eq. (|18p . is substituted into the power equation, 
Eq. p5|) . the following expression is obtained to first or- 
der in aiAT 



P 



1 A{l + aiAT) 
i?p,e {l + A{l + aiAT)y 



A 



i?prc (1 + A)2 



l + l-Aa,AT]u\ (29) 



This shows that an input voltage which is a pure 1. 
harmonic leads to dominant 0. and 2. harmonics in the 
power/temperature. Notice that if the thermistor resis- 
tance were temperature independent, that is ai =0, then 
Pa equals IP2I exactly but in general they are almost iden- 
tical, Po « 

It can further be seen that even if the voltage generator 
is purely single harmonic, the power can have a AuJc com- 
ponent from the interplay between the 2uj(. temperature 
variation and the 2cjo variation of the squared voltage. 
In principle further higher harmonics are generated this 
way. It can however be observed that if aiAT is small, 
or if A « 1 which can be chosen in the setup, the higher 
harmonics will be small (and decreasing with order). 

Eq. (P^l is not very practical for calculating the 
thermal-power terms needed for evaluating the thermal 
impedance from the temperature amplitudes obtained 
from Eq. since AT is needed. 

Alternatively the thermal power can be written as 

{U~V)V 



P = I{U -V)^ 



Rpi- 



(30) 



pre 



using Eq. (|14p . From this equation the power can be 
directly determined from measurable quantities. 

If only the Iwe component of the voltage is taken into 
account the following approximations are obtained 

P « ((C/i - Vi)Ei + c.c.){ViEi + c.c.) = 



4i?o 



-((C/i - V,)V{ + (C/i - V,)V,E2 + c.c.) (31) 
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leading to power coefficient given as 
1 



Po 



i{Ui-Vi)V{ + {U^ -V{)Vi) (32a) 



and 



Pa 



2-Rpri 



■(C/i - Vi)Vi. 



(32b) 



A more general expansion including higher harmonics is 
given in appendix [Bl 



The AC-method refined 



As mentioned earlier we have to do a refined analysis 
of the detected harmonics of the voltages. These extra 
terms stem from different sources: f ) The voltage source 
itself has higher harmonics and a DC-offset. 2) A sec- 
ond order term in the temperature dependence of the 
thermistor resistance. 3) The Joule-power compo- 
nent produced by the 2wo variation of the thermistor (as 
discussed in section IfV Ep . These effects are taken into 
account, as described below. 

Regarding point 1) we write the source voltage, U{t), 
more generally as 

1 



U4E4 + c.c.) 

(33) 

including up to the fourth order harmonics. A f/3 com- 
ponent from the source is to be expected since the source 
itself has an inner resistance. The Uq is a DC-offset that 
was present in this experiment but in principle should be 
eliminated. 

Regarding point 2) we expand the temperature depen- 
dence of the thermistor to second order as 



R = i?o(l + aiAT{t) + a2AT{tf). 



(34) 



Because R{T) follows Eq. (|19p over an extended temper- 
ature range a2 and ai are connected as 



ao = — 

9 T 



Oil 



(35) 



cryo 



Point 3) is taken into account by expanding V{t), P{t) 
and T{t) like U (t) in Eq. ([55]) up to the fourth harmonic 
components also. 

The full analysis can be found in the appendix [XJ includ- 
ing terms larger than 5 • 10~^C/i. Fortunately it leads to 
rather simple extensions of the Eqs. [28] 



^ V3iA + l)-iU:,+X3) 
aUi 



(36b) 



Xi and X2 depends on Tq, Ti, T2, T4, and b = [jfj) ^ - 
, and the expressions are given in Eq. (jA3P . 

The coupled Eqs. can be solved by iteration, us- 

ing a procedure equivalent to the one for the first order 
solution. 

We start by estimating A{T) from Eq. ()36ap . putting 
Xi = Tq = T2 ~ 0, and from this estimate ai and 012. 

A provisory thermal impedance Z2 is then found from Eq. 
P6bp . Ti = ZiPi and = Z4P4 can then be estimated 
since Zi(w) = 2'2(w/2) and ^4(0;) = Z2(2w), and To 
is found from the limit of Zq as described earlier. The 
perturbative terms, Xi^X^, can then be found and with 
these values in Eqs. (|36ap and (|36bp a better estimation 
of A (and hence ai and a2) and T2 is found. This process 
is iterated until convergence. 

The two important power components are calculated 
from the Vi's and Ui's as in Eq. ([32]) by 



Po 



and 



■J— (woVo + Irc {WiV{ + W2V2* + W3V3* + W4V4) 

Rprc V ^ 

(37) 



W0V2 + W2Va 



R 



pre 



- (WiVi + W^Vs + WiV4 + WsV^ + VKiKf) j (38) 

where W = U —V . The full set of power components are 
given in Eq. (|Bip . 



EXPERIMENTAL 



The measurements were performed using a custom built 
setup (Tgj . [20j . The temperature was controlled by a cryo- 
stat with temperature fluctuations smaller than 2mK (see 
Ref. [l^for details on the cryostat). The electrical signals 
were measured using a HP3458A multimeter in connec- 
tion with a custom-built frequency generator as sketched 
on Fig. [1] (see Ref. for details on the electrical setup, 
but notice that in this experiment we used the generator 
and multimeter in a slightly altered configuration). The 
electrical setup allows for measurements in the frequency 
(;/e) range from ImHz up to fOOHz. 

The used thermistor was of the type U23UD from 
Bowthorpe Thermometries, and was positioned in the 
middle of a sphere as shown on Fig. [2] The details of 
the sphere is as discussed earlier unimportant, and in 
this case the piezoelectric bulk transducer [l^ is used 
as such. The sphere is filled at room temperature with 
additional liquid in the reservoir. This additional liquid 
ensures that the sphere stays completely filled as it is 
cooled to the relevant temperature. 
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Cryostat 



FIG. 2: Schematic illustration of the measuring setup. The 
thermistor acting as combined heater and thermometer is 
placed in a sphere filled with the liquid. The sphere is it- 
self placed in a cryostat. 



The liquid investigated was a five-ring polyphenyl ether 
(5-polyphenyl-4-ether) the Santovac@ 5 vacuum pump 
fluid (CAS number 2455-71-2). The hquid was used as 
received, without further purification. 



Data were taken as frequency scans at constant temper- 
ature, with the liquid in thermal equilibrium. The direct 
output from the frequency generator was measured under 
electrical load, at each sample temperature. This elimi- 
nates any long term drifts in the frequency generator out- 
put voltage. The measurements were furthermore carried 
out at two different input amplitudes, {\Ui\) « 4.9V and 
{\Ui\) w 2.9V, allowing for a test of the inversion algo- 
rithm and linearity of the final thermal response. If not 
stated otherwise all results presented in this paper are 
from the larger of the two amplitudes. To test for re- 
producibility and equilibrium some of the measurements 
were retaken during the reheating of the sample. The 
outcome of the measurements are the complex ampli- 
tudes Ui{vc) and Vi(i'c) for the harmonics of the voltages 
U and V. 



Table U gives the key experimental parameters for the 
measurements at the high input amplitude. 



VI. DATA ANALYSIS 

In the following we denote the measured thermal 
impedance, Z2^ by simply Z dropping the subscript 2. 
The subscripts were practical in the technical descrip- 
tion of the detection method indicating that the thermal 
impedance is found at a thermal frequency, u, which is 
the double of the electrical frequency, Wg. But in dis- 
cussing the thermal properties of liquids only the ther- 
mal frequency lo = 2t:v is relevant. The temperature 
we refer to at which a measurement is done is the mean 
temperature at the bead surface, T = Tq -\- Tcryo- 

The thermal impedance has been measured at two differ- 
ent input voltage amplitudes, Ui, of the voltage divider 
of Fig.[T]and thus at two different power amplitudes. The 
liquid thermal response, the temperature, is expected to 
be linear in the power amplitudes and thus the thermal 
impedance should be independent of the power ampli- 
tude. 

The thermal impedance has been calculated from the 
measured voltage harmonics using both the first-order 
method described in section IIVDI and using the higher- 
order method described in section HV Fl 

The top part of Fig. [3] shows the frequency depen- 
dent relative difference between the calculated thermal 
impedance from the two amplitudes. The results from 
the two methods (the first-order and higher-order meth- 
ods) for three selected temperatures are shown in the 
figure. In the lower part the average relative difference 
is shown as function of temperature. It is seen that the 
higher order method generally reduces the difference sig- 
nificantly. Only in a very small temperature interval does 
the two methods give equally good results. This demon- 
strates the necessity of taking into account the higher- 
order terms of the elaborated analysis (as described in 
section [IVF| . 

Fig. m shows the real and imaginary part of the thermal 
impedance, Z, in a log-log plot at 295. 6K. The mea- 
sured impedance has more features than predicted by the 
simple expression for the thermal impedance of spherical 
effusion Eq. Q. For example the high frequency behav- 
ior should be characterized by a line of slope 1/2 but 
this is seen not to be the case. The additional features 
come from thermal properties of the thermistor bead, 
the influence of which we are going to study below. We 
stress that although Z shows dispersion and has a real 
and imaginary part this has nothing to do with liquid 
relaxation — which is absent at this temperature — but 
is only a consequence of heat diffusion in spherical ge- 
ometry. However relaxation should of course affect the 
thermal impedance. This is shown in Fig. [5] where Z at 
252. 7K and 256. 7K are compared. The glass transition 
is seen to give only a slight perturbation of the shape 
of the Z-curve. In order to obtain any reliable specific 
heat data from the thermal impedance it is thus impera- 
tive to have an accurate model of the thermal interaction 
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FIG. 3: Comparison of the thermal impedance, Z, obtained 
from the measurements at two different ampUtudes ((C/i) = 
4.9V and 2.9V), and from using the first-order model as given 
in Eq. (|28p and higher-order model as given in Eq. (|36|) (as 
described in section HV Dl and IIV F[) . 

(A) the frequency dependent relative difference between the 

two methods (|Zamp=2.9V - .^amp=4.9v|/|^amp = 2.9v|). Full 

lines: the difference for the higher-order model. Dashed 
lines: the difference for the first-order model. Cryostat tem- 
peratures are 293. 5K (blue), 261K (green) and 250. 5K (red). 

(B) shows the frequency averaged relative difference 

({|^amp = 2.9V — Z, 



FIG. 4: Comparison of measured thermal impedance, Z, with 
the result from fitting the full model of liquid and thermistor 
(Eqs. (O, ((42]) and (gS])), and the simple model of Znq (Eq. 
(|C2[l 'l. The parameters used in evaluation of the simple model 
are the ones found from the full analysis. Data are taken at 
room temperature (T — 295. 6K). Dots: measured points. 
Full lines: full model of liquid and thermistor bead. Dashed 
lines simple model including only the liquid. 



between liquid and thermistor bead. 

The bead (see Fig.lH]) consists of a core of radius tq, where 
the temperature amplitude 5T is actually measured and 
the heat current P generated. This core is surrounded 
by a glass capsule of outer radius ri, at which the bead 
is in contact with the liquid. The heat flow P{ri) out 
through the surface at ri is in general different from the 
heat flow P{ro) out through the surface at tq. Also the 
temperature ST{ri) can be different from (5T(ro). Now 
the liquid thermal impedance is 



Zliq — 



Pin) 



(39) 



whereas the measured thermal impedance is rather 

ST , , 

Z=—. (40) 

The capsule layer has a heat conductivity Af, and a spe- 
cific heat Cft. The bead is a solid for which in general 
the difference between the isobaric and isochoric specific 
heats are small and they are frequency independent. The 
heat diffusion through this shell is thus well described [l^l 



4.9v|)./(|^amp=2.9v|>.). Squares the ^ ^ transfermatrix, T"^ = T"^(Afc, c^, n, tq), 



difference for the higher-order model. Circles the difference 
for the first-order model 



<5T(ri) 
P(ri)/(i^) 



-ith 



(41) 
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This implies that the thermal impedance Znq = Zr-^ at 
radius ri is transformed to a thermal impedance Zj-g at 
radius tq via 



Zr, 



^12 



T^^iujZ^ 



liq 



(42) 



A small part of generated power is stored in the core of 
the bead. The heat capacity Ccoro of the core is taken 
to be 4/37rro'^Cb. The measured thermal impedance Z is 
thus related to Zro by 



1 _ 1 

— — IWOcorc + 

Zj Zjrr 



(43) 



The model relating Zuq to the measured Z described by 
Eqs. ([9]), (|42|) and (|43l) are depicted in the electric equiv- 
alent diagram in Fig [SI 

In order to get from the measured Z one has to 
know the parameters Tq , ri , A;, and characterizing the 
thermal structure of the intervening thermistor bead. 

This can in principle be done by fitting Z of the model 
obtained by combining Eqs. ©, (gl]), and (|33]) to the 
measured frequency dependent thermal impedance. An 
electric equivalent network of the model is depicted in 
Fig. [S] This model has six frequency independent intrin- 
sic parameters, namely tq, ri, A;,, Cf,, A and c;, if the fit is 
limited to the temperature region where c; can be consid- 
ered frequency independent. By careful inspection of the 
equations it can be found that the number of indepen- 
dent parameters that can be deduced from a fit is only 
five. The five parameters can be chosen in a number of 
ways, of which we have chosen 




-1 1 

log^o(v) [Hz] 



FIG. 5: Comparison of thermal impedance, Z, at two temper- 
atures, showing the small difference due to the frequency de- 
pendent specific heat of the liquid, ci(y). Lines: T — 256. 9K. 
Dots T = 252. 8K. 



Glass capsule 
Heat generator \^ 



2^; 



Xb ri Cb 47rAri 



\ 



(44) 



This set of variables has the advantage that it minimizes 
the number of parameters which mix bead and liquid 
properties. Here c is the only one such mixed parameter. 
The expressions for the transfermatrix and Ziiq in terms 
of these parameters is given in appendix [Cl 

The five parameter model was fitted to the measured 
thermal impedance, Z, by using a least square algorithm. 
The results are shown on Fig. [7| and [8l Measurements 
were taken both going downwards to the glass transition 
and upwards again to room temperature. In the glass 
transition range where the liquid specific heat, c;, is ex- 
pected to be complex and frequency dependent these fits 
of course do not work. This is the reason for the scatter 
of data at low temperature. This difference between the 
high and low temperature regimes can also be seen di- 
rectly in the error from the fit as illustrated on the lower 
part of Fig. [51 

The solid line in Fig. [3| shows the result of such a fit to Z- 
data at room temperature (295. 6K). The model captures 




Thermistor 



Liquid 




-hq 



I'd, /"i, Cf), A5 



Cl, A 



FIG. 6: Schematic thermal network model of the thermis- 
tor bead and liquid. The thermal current P and tempera- 
ture ST axe indicated at the inner radius, ro, at the border 
between the heater and the glass capsule and at the outer 
radius, ri, at the border between the glass capsule and the 
liquid. The thermal impedance of the liquid is defined as 
Zuq ~ 5T{r\)/ P{ri). The thermal impedance seen from the 
thermistor is Z = 5T{ro)/P. The two are connected by Eq. 
(|42p and (|43p. utilizing the thermal transfermatrix Eq. (|C3[l . 
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FIG. 7: Results from fitting the measured thermal impedance, 
Z, to the 5-parameter thermal network model (Fig. [S| of the 
spherical heat effusion experiment Eqs. ((9}, (|42|) and (|43|) . 
This figure shows (a) r = ro/ri , (b) c = cj/cft and (c) Tt, 
(the remaining parameters are shown on figure [8]). 
Open squares (red): Data taken going down in tempera- 
ture. Open circles (green): Data taken going up in tem- 
perature. Closed circles (blue): Low temperature data 
including points take going down and up. Line: Interpo- 
lated temperature dependencies of f , c and Tb found based on 
data from above the liquid glass transition temperature range, 
used for extrapolating down in the region where the liquid re- 
laxation sets in (the functional form and fitting results are 
shown). 
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FIG. 8: Results from fitting the measured thermal impedance, 
Z, to the 5-parameter thermal network model (Fig. [6]) of the 
spherical heat effusion experiment Eqs. Q, (|42p and (|43p . 
This figure shows (a) rj, (b) Zuq,o, and (c) the relative error 

from the fit, ar = '\J ^Yl^\z^"'d'lf^^ ^^^^ remaining parame- 
teres are shown on figure [7| 

Open squares (red): Data taken going down in tempera- 
ture. Open circles (green): Data taken going up in tem- 
perature. Closed circles (blue): Low temperature data 
including points take going down and up. Line: Interpo- 
lated temperature dependencies of r;, and Zuq.o found based 
on data from above the liquid glass transition temperature 
range, used for extrapolating down in the region where the 
liquid relaxation sets in (the functional form and fitting re- 
sults are shown). 



the features in the frequency dependence of the measured 
thermal impedance very precisely. Also in Fig.[4]is shown 
the thermal impedance as it would have been if the prop- 
erties of the thermistor bead itself had been negligible, 
i.e. if Z was identical to Zn^ of Eq. ([9]) We note that 
especially the high frequency behaviors are different but 
also the limiting value at low frequency is higher for the 
measured Z than for Z^q. This is due to the additional 
thermal resistance in the bead capsule layer. 

In table [ill we report the results for Tb and c obtained 
at room temperature. Four values are given for each 
quantity, they are taken at the two input amplitudes and 
from an initial measurement at room temperature and a 
measurement taken at the end of the temperature scan. 
The method is seen to give very reproducible results. 



The five parameters are extrapolated down to the regime 
of frequency dependent c;, by the functions indicated on 
the figures, r is considered being temperature indepen- 
dent, c, r and are to a good approximation linear 
in temperature over the investigated temperature range. 
Ziiq is found to be very well fitted by a second degree 
polynomial. 

The extrapolated values of the five parameters can then 
be used for calculating Zuq from the measured thermal 
impedance, Z inverting the relations Eqs. (|43|l and (|42| . 
In Fig. [9] the thermal impedance of the liquid, Zuq, and 
the raw measured thermal impedance, Z, are shown for 
data taken at room temperature (= 295. 6K). 

As discussed earlier effusion from a plane plate can only 
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TABLE II: Heat conductivity and specific heat for the hquid (A and c;), as defined in relation to Fig. 1101 and |12I Fitting 
parameters characterizing the bead (rt and c) from the fits to the five-parameter mode. The results are at room temperature 
(295. 6K), and taken using two input amplitudes |l7i|. The two scans are the initial scan and a scan taken after the temperature 
scan coming back to room temperature. The results indicates the reproducibility of the method. 
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FIG. 9: Comparison of the raw measured thermal impedance, 
Z, and final liquid thermal impedance, Zuq, corrected for the 
influence of the thermistor encapsulation. Data are taken at 
room temperature (T = 295. 6K), and are equivalent to the 
data shown on Fig. |4l The liquid thermal impedance is trun- 
cated at high frequencies {u > 10^'^ Hz) as the bead becomes 
dominant at higher frequencies, leading to unphysical behav- 
ior of the calculated liquid impedance. Full lines: Raw ther- 
mal impedance Z. Dashed lines: Liquid thermal impedance 

Zuq. 



give the efFusivity e = VXc. Thus the specific heat cannot 
be found absolutely. In contrast effusion from a sphere is 
able to give both the heat conductivity and the longitudi- 
nal specific heat absolutely even if the latter is frequency 
dependent. However, the heat conductivity has to be 
frequency independent to allow for this separation. The 
separation is possible because Z\iq as given in Eq. ([9]) has 
the low frequency limiting value Znq Q = l/(47rAri). The 
convergence is rather slow involving the squareroot of the 
frequency. However, if we look at the reciprocal quantity. 



the thermal admittance Ynq = l/Znq 

Hiq = 4TT\n (l + v/«wri2Q/A) (45) 

we observe that when ci is frequency independent the 
squareroot frequency term can be cancelled by taking 
the difference between Y^[^ and Y^['^ 

y - Y" 

A = ^. (46) 

This is valid even in the glass transition range at suffi- 
ciently low frequencies where ci reaches its non-complex 
equilibrium value. Finally, using the found values for A 
the frequency dependent longitudinal specific heat is cal- 
culated from 

Afl^-l)' (47) 
lujri \47rAri / 

However the structure of the five parameter model only 
allows for determining f but not rg and ri independently 
from the fit to data. To obtain the absolute value of A and 
ci one of the radii are needed. Measured values of ri are 
given in table |T1 As indicated the bead is not completely 
spherical, and it is therefore needed to determine which 
value of ri to use. 

According to the specifications for Santovac@ 5 given 
by Scientific instrument service, Inc (2l| the value of A 
at 20°C is 0.1330W/(K-m). They likewise report the 
isobaric specific heat pr. mass Cp/ p = 0.35cal/(g • K) and 
density p = 1.204g/cm'^ at 20° C giving a specific heat pr. 
volume of Cp = 1.76 • 10^J/(K • m'^). At these tempera- 
tures the isobaric and longitudinal specific heats are iden- 
tical, we have therefore chosen to use a value for ri which 
matches q, within the errors, to the reference value. This 
"optimized" value is very close to the ri^Moasurcd, long edge 
value given in tableU This choice gives a value of A which 
is 3.5% too big compared to the literature value. 

Fig. [To] shows how the temperature dependence of the 
heat conductivity, A, of the sample liquid is found. In 
Fig. [11] is shown the final frequency dependent specific 
heat for a number of temperatures. Table iHl reports the 
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FIG. 10: Thermal conductivity of the liquid, (a): as 
function of frequency, plotted for a selection of temperatures 
(same temperatures as on Fig. llll with colors indicating same 
datasets). According to Eq.|35]the thermal conductivity is the 
low frequency limit of this quantity. The points on the lines 
indicate the data points used for defining the low frequency 
limit (taken as an average over these points), (b): Thermal 
conductivity of the liquid defined from the average over the 
y'- —y'-' 

low frequency limit of ''^^^ . Squares corresponds to the 
curves shown on the upper figure, and dots to the remaining 
datasets. 
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FIG. 11: Real and imaginary part of the longitudinal specific 
heat of the liquid, cj, at a number of the investigated temper- 
atures. The frequency axis is truncated at O.SHz as the signal 
above this frequency is dominated by the inner-structure of 
the thermistor. Temperatures are 295.6K, 271. IK, 261. OK, 
256.9K, 254.9K, 253.8K, and 252.8K. 



glass transition. At 252. 7K q decreases from Cq = 
ci{uj 0) = 1.75-106j/(K-m3) to = ci{uj oo) ^ 
1.30 • 10^ J/ (K • m"^). The temperature dependency of cq 
above the glass transition range is shown in Fig. [TH The 
temperature dependency of the characteristic relaxation 
time, T, is shown in Fig. [131 The high temperature data 
are seen to show spurious frequency dependency in both 
the real and imaginary part of q. This may be caused 
by small systematic errors in the data processing or data 
acquisition that we have not been able to trace yet. It 
may also be due to the small eccentricity of the bead. 
Due to the functional relationship, Eq. ((21) between Znq 
and ci a small relative error \dZiiq/ Znql in the thermal 
impedance may propagate to a large relative error \dci jci \ 
in the specific heat. As a function of the dimensionless 
Laplace-frequency s = iujTi one finds 



room temperature values for c; and A for two tempera- 
ture amplitudes, and data reproduced after a tempera- 
ture scan down and up again. The values indicate the 
reproducibility of the method. 



VII. DISCUSSION 

The derived longitudinal specific heat displayed in Fig. 
111! shows the expected relaxation phenomena at the 



dci , dZxiq _ q I + 7g (-43^ 

Cl Zliq ^/S 

This factor is four at the characteristic heat diffusion fre- 
quency (s = 1) increasing to a factor of 22 at 0.01 times 
this frequency. Thus the specific heat is difficult to get 
reliable more than 2 decades below the characteristic dif- 
fusion frequency. On the other hand one decade above 
this frequency the measured Z is dominated by the ther- 
mal structure of the thermistor bead that can only to 
some extent be reliable modelled . 
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FIG. 12: DC specific heat, co, of the iiquid for temperatures 
above the liquid glass transition temperature range. The spe- 
cific heat is found as the average over the real part of ci in 
the frequency range 10~^Hz-lHz. cq corresponds to the DC 
value of Cp as can be seen from Eq. ([3| and Eq. Open 
squares (green): Data taken going down in temperature. 
Open circles (red): Data taken going up in temperature. 



VIII. CONCLUSION 



From a theoretical point of view determination of the 
frequency dependent specific heat by a thermal effusion 
method based on a sphere has many ideal features. The 
surface of a sphere — being closed — has no boundary 
and thus no associated boundary effects like a finite plane 
plate. The thermomechanical coupling problem can be 
treated analytically i.e. the influence of the increas- 
ing dynamic shear modulus and the mechanical bound- 
ary conditions are adequately taken into account. Al- 
though the thermomechanical problem can also be solved 
in the onedimensional unilateral case it is not so obvious 
whether it really applies to the experimental finite plate 
realizations. It is of course preferable if thermal experi- 
ments can give two independent thermal properties e.g. 
the specific heat and the heat conductivity or the effu- 
sivity and diffusivity. Effusion experiments in different 
geometries give the effusivity only unless a characteris- 
tic length scale comes into play with the heat diffusion 
length. For the infinite plane there is no such length 
scale whereas in the spherical case this length scale is the 
radius of the heat producing sphere. This advantage is 
however also the weakness of the method since it limits 
the practical frequency window to be studied. 
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FIG. 13: Characteristic relaxation time, r, for the longitudi- 
nal specific heat, ci (oj). The time if found from a fit to 2. order 
polynomial to a few point around the maximum of the imag- 
inary part of Ci and corresponds to r = l/(27ri^ip), where i/ip 
is the frequency of maximum loss. Open squares (green): 
Data taken going down in temperature. Open circles (red): 
Data taken going up in temperature. 
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Appendix A: The AC-method to higher order 



In this appendix we analyze which higher-order terms must be included in order to get a relative accuracy of 1 • 10 ^ 
on the measured temperature amplitude from the Fourier components of the measured voltages of the voltage divider 
(see Fig. [T|). 

In the following we will include up to fourth harmonics in all the relevant quantities, hence we write 

(Ala) 
(Alb) 
(Ale) 
(Aid) 



V{t) = 


^{Vo + ViEi+V2E2^ 


- V3E3 - 


f ViE4 4 


- c.c.) 


U{t) = 


i [Uo + UiEi + U2E2 


+ U3E3 


+ UiEi 


+ c.c.) 


Pit) - 


]^{Pq + PiEi + P2E2 - 




+ P4E4 - 


fee.) 


AT{t) = 


i (To + T^Ei + T2E2 4 


- T3E3 - 


f T4E4 4- 


- c.c.) . 



Combining the expansion of the thermistor resistance to second order, Eq. 
Eq. (fT6|) . the voltage across the preresistor becomes, to second order in AT, 

1 



V(t) = (1 - aAT + bAT^) U 



with the voltage divider equation, 

(A2) 



with A = Ra/Rpyc, 



^and6: 



1+A 



Aa2 
1+A- 



If this expression is explicitly calculated using the above expansions of U and T (Eq. (jAlbp and (jAld|l ) a large number 
of terms is obtained. In the following we make a numerical inspection of which terms are most significant, in a worst 
case, to reduce the number of terms included. 

In order to calculate the order of magnitude of the terms, an estimate of the size of the components of U and AT 
are needed. In Fig. [T3] we show the relative magnitudes of the different components of the input voltage, U, of the 
voltage divider to C/i. In table Hill we summarize the upper limits of these ratios. In Fig. [15] we show the relative 
magnitude of the different components of the power calculated from the measured Ui and Vi components. In table HVl 
we summarize the upper limits on these ratios. We estimate the relative size of the temperature amplitudes as equal 
to the relative sizes of the power amplitudes. In table |V] we give characteristic values of the quantities describing the 
thermistor and the resistors. 

Combing the estimates given in table |Vl with Eq. ((26l) . we observe that a change in r2 of 1 • 10^'*K corresponds to a 
change in V3/U1 of 5 • 10~^. Hence only terms in the expansion of Eq. (|A2[) which are larger than 5 • 10~^[/i will be 
included. 



By explicit substitution of Eq. (jAlbp and (jAldp into Eq. (jA2| and insertion of the estimates from table lllli IIVI and 
rVl the following expressions are found for the first and third harmonics on the voltage, V, when disregarding terms 
smaller than 5 • 10~^J7i 



Vi = 



1 



A + l 



Ui - aToUi - -aT2U* + Xi 



(A3a) 



and 



V3 = 



A+l 



--aT2Ui + U3 



X3 



(A3b) 



with 



Xi = ^aTiUn + bT^Ui + -bT2T;Ui + bToT2U* 
X3 = -^aTiU* + hT^U* + bToT2Ui+bToTiU*. 



(A3c) 
(A3d) 
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Appendix B: Power terms to higher order 

The relevant set of power components are calculated from Eq. ([50)1 as 

{u-v)v 

^prc 

using the expansions in Eq. ([Hal and ((XIb| to be 

Po = (woVo + ^Re {WiV{ + W2V^ + W^V^ + ^^4^4*)) (Bla) 

Pi = (WiVq + ViWa + \ {W2V1 + 1^1*^2 + M^3^2* + M^2*^3 + W3*TAj + 1^4 K,*)) (Bib) 

P2 = (woV2 + W2V0 + I {WiVi + W^Vs + W;Vi + 1^3 n* + WiV^)] (BIc) 

-Kpro V ^ / 

P3 = (woV^ + W^Vq + i {W1V2 + M^2^i + W^Vi + WaV*)^ (Bid) 

Pi = (Wo 14 + WiVo + \ {WiVz + W2V2 + WzVi)] (Ble) 

where W = U ~ V . 



Appendix C: The "Ave parameter" formulation of the thermal transfer model 



Eq. (I9| , (|41l) and (|42|) are in the following written in terms of the five variables 

2^1 2Cb . ro ^ ci 1 

^i=^iT' '^b^r^—, r = —, c=— , Ziiq,o = — : — . (CI) 

A Ab ri Cb AnXri 

The thermal impedance, Eq. ([9|), is given as 

= (TTT^ ^^^^ 

The components of the thermal transfer matrix [l3| in Eq. (j4T|) are given as 



Tl"^ = f cosh (VioJnil - f)) + i sinh (\/iZ^(l - f)) (C3a) 



Ti*2h = -Zuq.o-^^sinh(Vi^(l-r)) (C3b) 
Ti r 



rpth _ 1 '^i 1 

^21 



(1 — iujTbr^ sinh (•\/zar7^(l — f)) — \/*^''"6(l ^ f) cosh (•\/wIr7^(l — f)) (C3c) 



Tn'o = T cosh (^/IZJ^(l - f)) 7^sinh(Vi^(l - r)) (C3d) 

r ry/iujTb 
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\09,,{v^) [Hz] 



TABLE III: Maximum of relative magnitude of the compo- 
nents of input voltage U to the dominant first harmonic Ui, 
as shown on Fig. [14] (for u < 40Hz). 

~Uo/Ui U2/U1 U3/U1 Ua/Ui^ 
2-10'^ IQ-^ 10'^ 3 ■10''"' 

°'Ua/U\ was estimated from an additional measurement where U4, 
was included. 



FIG. 14: Relative magnitude of the components of input volt- 
age U to the dominant U\ component. (|?7i|) = 4.909V. See 
also table Hill 

TABLE IV: Maximum of relative magnitude of the compo- 
nents of the power as shown on Fig. 1151 (for v < 40Hz). 

P0/P2 P1/P2 P3/P2 P4/P2 " 
1 10"^ 10"^ 10"^ 




FIG. 15: Magnitude of the components of the power to the 
dominant P2 component. {jf/i|> = 3.927V, Tcryo = 249.5K. 
See also table HVl 
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